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Abstract. We introduce a random walk in random environment associated to an under- 
lying directed polymer model in 1+1 dimensions. This walk is the positive temperature 
counterpart of the competition interface of percolation and arises as the limit of quenched 
polymer measures. We prove this limit for the exactly solvable log-gamma polymer, as 
a consequence of almost sure limits of ratios of partition functions. These limits of ra- 
tios give the Busemann functions of the log-gamma polymer, and furnish correctors that 
solve a variational formula for the limiting free energy. Limits of ratios of point-to-point 
and point-to-line partition functions manifest a duality between tilt and velocity that 
comes from quenched large deviations under polymer measures. In the log-gamma case 
we identify a family of ergodic invariant distributions for the random walk in random 
environment. 
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1. Introduction 

In directed polymer models the definition of weak disorder is that normalized point-to- 
line partition functions converge to a strictly positive random variable. In strong disorder 
these normalized partition functions converge to zero. Weak disorder takes place only in 
dimensions 3+1 and higher and under high enough temperature, while lower dimensions 
are in strong disorder throughout the temperature range. (See [3, 5, 6, 7, 8, 12, 20] for 
reviews and some key results.) 

We work in 1+1 dimensions with the explicitly solvable log-gamma polymer. We show 
that ratios of both point-to-point partition functions and tilted point-to-line partition func- 
tions converge almost surely to gamma-distributed limits. Out of this basic fact we derive 
several consequences. 

(i) Limits of ratios of partition functions give us limits of quenched polymer measures, 
both point-to-point and point-to-line, as the path length tends to infinity. The limit pro- 
cesses can be regarded as infinitely long polymers. Technically they are random walks in 
correlated random environments (RWRE). When we average over the environment, this 
RWRE has fluctuation exponent 2/3, in accordance with 1+1 dimensional Kardar-Parisi- 
Zhang universality. This polymer RWRE is also a positive temperature counterpart of a 
competition interface in a percolation model. (This terminology comes from the idea that 
percolation models are zero-temperature polymers. Remark 2.1 below explains.) For the 
RWRE we identify a family of stationary and ergodic distributions for the environment as 
seen from the particle. The averaged stationary RWRE is a standard random walk. 

(ii) Logarithms of the limiting point-to-point ratios give us an analogue of Busemann 
functions in the positive temperature setting. Busemann functions have emerged as a 
central object in the study of geodesies and invariant distributions of percolation models 
and related interacting particle systems [1, 4, 11, 16, 21]. Our paper introduces this notion 
in the positive temperature setting. We show how Busemann functions solve a variational 
problem that characterizes the limiting free energy density of the log-gamma polymer. 

A theme that appears more than once is a familiar large deviations duality between the 
asymptotic velocity of the path under polymer distributions and a tilt introduced into the 
partition function and probability distribution. In this duality the mapping from velocity 
to tilt is given by the expectation of the Busemann function. In particular, this duality 
determines how limits of ratios of point-to-point and tilted point-to-line partition functions 
match up with each other. 

A word of explanation about our focus on the log-gamma polymer. The ultimate goal 
is of course to find results valid for a wide class of polymer models. We could formulate 
at least some of our results more generally. But the statements would be complicated and 
need hypotheses that we can presently verify only for the log-gamma model anyway. For 
general polymers, just as for general percolation models, we cannot currently prove even 
mild regularity properties for the limiting free energy. Thus we chose to focus exclusively 
on the log-gamma model (except for the general discussions in Sections 2 and 5). 

We expect that much of this picture can eventually be verified for general 1+1 dimen- 
sional directed polymers. Our hope is that this paper would inspire such further work. For 
example, it is clear that the solution of the variational formula for the free energy in terms of 
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Busemann functions works completely generally, once a sufficiently strong existence state- 
ment for Busemann functions is proved. Busemann functions with tractable distributions 
are an essential feature of the exact solvability of the log-gamma polymer. They can be 
used to construct a shift-invariant version of the polymer model, which was earlier used for 
deriving fluctuation exponents and large deviation rate functions [15, 27]. 

The log-gamma polymer is a canonical model in the Kardar-Parisi-Zhang universality 
class, in the same vein as the asymmetric simple exclusion process and the corner growth 
model with geometric or exponential weights [9, 17, 23, 28, 29]. It was introduced in [27] 
and subsequently linked with integrable systems and interesting combinatorics [2, 10, 22]. 
These exactly solvable models are believed to be representative of what should be true more 
generally. 

Organization of the paper. The paper is essentially self-contained. One exception is that 
in Section 5 we cite variational formulas for the free energy from [24, 25]. Here is an outline 
of the paper. 

Section 2. Introduction of the polymer RWRE in a general context as the positive 
temperature counterpart of the competition interface of last-passage percolation. 

Section 3. Introduction of the log-gamma polymer. The shift-invariant log-gamma poly- 
mer is formalized in the definition of a gamma system of weights. 

Section 4. Limits of ratios of point-to-point partition functions for the log-gamma poly- 
mer. 

Section 5. Busemann functions are constructed from limits of ratios of point-to-point 
partition functions and used to solve a variational formula for the limiting free energy. 
Duality between tilt and velocity. 

Section 6. Limits of ratios of tilted point-to-line partition functions for the log-gamma 
polymer. Duality between tilt and velocity appears again. 

Section 7. Limits of ratios of partition functions yield convergence of polymer measures 
to the polymer RWRE. The limit RWRE has fluctuations of size n^/^ under the averaged 
measure. 

Section 8. A stationary, ergodic distribution for the log-gamma polymer RWRE. 
Section 9. Several auxiliary results, including a large deviation bound for the log-gamma 
polymer and a general ergodic theorem for correctors. 

Notation and conventions. N = {1,2,3,...} and = {0, 1,2,...}. For n E N, [n] = 
{1,2,..., n}. On the £^ norm is |x|i = + |x2|, and inequalities are coordinatewise: 
(xi,X2) < (2/1,2/2) if < Vr for r G {1,2}. Standard basis vectors are ei = (1,0) and 
62 = (0,1). Our random walks live in and admissible paths x, = (xfc)^^Q have steps 
Zk = Xk — Xk-1 € {61,62}. Points of are written as u,v,x,y but also as {m,n) or 
Weights indexed by a single point do not have the parentheses: if x = then r/^. = r]ij. 

For u < V in Z^, 11^ ,^ is the set of admissible paths from xq = u to = v. Limit 

velocities of these walks lie in the simplex U = {(n, 1 — u) : u & [0,1]}, whose (relative) 
interior is denoted by int U. Shift maps T-n act on suitably indexed configurations w — {wx) 
by {T^w)x = Wv+x- IE and P refer to the random weights (the environment), and otherwise 
E^^ denotes expectation under probability measure /i. The usual gamma function for /o > 
is r(p) = xP~^e~^ dx, and the Gamma(/>) distribution on is T{p)~^xP~^e~^ dx. 
\E'o = r'/r and \I'i = \I'q are the digamma and trigamma functions. 
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The reader should be warned that several different partition functions appear in this 
paper. They are all denoted by Z and sometimes with additional notation such as Z. It 
should be clear from the context which Z is meant. Each Z is a sum of weights W{x}j 
of paths X, from a collection of nearest-neighbor lattice paths. Associated to each Z is a 
polymer probability measure Q on paths: Q{x.} = Z-'^W{x.). 



2. The polymer random walk in random environment 

In this section wc introduce a natural random walk in random environment (RWRE) 
associated to an underlying directed polymer model. This walk appears when we look for 
the positive temperature counterparts of the notions of geodesies and competition interface 
that appear in last-passage percolation. Percolation and polymers arc discussed in this 
section in terms of positive weights, without specifying probability distributions. 

2.1. Geodesies and competition interface in last-passage percolation. We give a 

quick definition of last-passage percolation, also known as the zero-temperature polymer. 
Let {ujx '■ X G Z^} be a collection of positive weights. For u < v m. let n^^^ denote 
the set of admissible lattice paths x, = (xj)o<i<n with n = |t; — n|i that satisfy xq = u, 
Xi Xi~\ e {61,62}, = V. The last-passage times are defined by 

\v-u\i 

(2.1) Gu,v = max ^ w^;., u<v in 

1=1 

A finite path (xi)o<i<n m ^u,v is a geodesic between u and v if it is the maximizing path 
that realizes G«,^, namely, Gu,v = Y17=i ^xi- Every subpath of a geodesic is also a geodesic. 
Let us assume that no two paths of any length have equal sum of weights so that maximizing 
paths are unique. This would almost surely be the case for example if the weights are i.i.d. 
with a continuous distribution. 

It is convenient to construct the geodesic from n to v backwards, utilizing the iteration 

(2-2) Gu^x ~ Gii^x—ei V Gu^x—e2 ~^ ^x- 

Start the construction with x„ = v. Suppose the segment {xk, Xk+i, ■ ■ ■ , Xn) of the geodesic 
has been constructed. If Xk > u coordinatewise, set 



(2.3) Xk-i 




if Gu^xk—ei G^UjXfc— 62 5 
if Gu,xk—ei < Gu,xk-e2- 



If Cr for either r = 1 or r = 2, then define the remaining segment as {xq, . . . ,Xk) = 

(u + ie3-r)o<i<k- 

For a fixed initial point u € Z^_, the geodesic spanning tree Tu of the lattice u + Z^ is 
the union of all the geodesies from u to v, v e u + Z'^. 

The competition interface (p = {(pk)kez+ is a lattice path on Z^ defined as a function of 
{Go,v}vez'^- starts at (^0 = and then chooses, at each step, the minimal G- value: 

(^2 4) ^ \vk + ei if Gc^^+ei < Go,^^+e2, 

\'Pk + e2 if Go,<^j,+ei > Go,¥'fc-|-e2- 
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The relationship between To and ip is that 99 separates the two subtrees 7o,ei , 70,62 of To 
rooted at ei and 62- Since every lattice path from has to go through either ei or 62, 
To = {0} U To,ei U To,e2 as a disjoint union. For each n G Z_(_, is the unique point such 
that l^^nli = n and for r G {1,2}, + : A: G N} C To,er- Note that we cannot say 
which tree contains ipn, unless we know that ipn — V^n-i = &r in which case G To,er- If 
we shift if by (1/2, 1/2) then it threads exactly between the two trees (Figure 1). 

The term competition interface comes from the interpretation that 7o,ei and To, 62 are 
two competing clusters or infections on the lattice [13, 14]. The model can be defined 
dynamically. The clusters at time t G M+ are To^erii) = {v € To,6r '■ Gq^v < 
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FlGURE 1. The competition interface shifted by (1/2,1/2) (solid line) separating 
the subtrees of To rooted at ei and 62- 



2.2. Geodesies and competition interface in positive temperature. We switch to 
a positive temperature polymer. Let {Vx}x^i?^ be positive weights. Define point-to-point 
polymer partition functions for u <v vl\ Z^ by 

1)— n| 1 

(2.5) z,,. = E n ^-t' 

and the polymer measure on the set of paths n„ by 

^ \v-u\-i 

(2.6) QuA^.) = ^ n ^-T'' ^. en,,.. 

' 1 = 1 

Our convention is to use reciprocals Y^^ of the weights in the definitions. The reason is 
that this way the weights in the log-gamma polymer are gamma distributed and features 
of the beta-gamma algebra arise naturally. 
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Remark 2.1. A conventional way of defining polymer partition functions is 
(2.7) e'^^l-'^-- 

with an inverse temperature parameter < /3 < oo. In the zero-temperature limit 
f3~^ log Zu,v — > Gu,v as /3 — > oo, and the polymer measure Qu,v concentrates on the ge- 
odesic (s) from u to V. This is the sense in which last-passage percolation is the zero- 
temperature polymer. See Remark 3.2 below for this point for the log-gamma polymer. 



We implement noisy versions of rules (2.3) and (2.4) to define positive temperature 
counterparts of geodesies and the competition interface. 

Fix a base point u £ and define a backward Markov transition kernel V" on the 
lattice u + by tr"(u, u) = 1, and 

(2.8) V"(x,x-e^) = y - = ^ —f , for rG {1,2}, 

if both X and x — lie in u + Z^. The middle formula above gives the correct values on the 
boundaries of n + Z^ where there is only one admissible backward step: '^"(n + ier,u + 
{i - l)er) = 1 for i > 1 and r € {1, 2}. 

For a path x, S ^u,v comparison of (2.6) and (2.8) shows 



\v—u\ 1 

TT [Xi , X-i 
1=1 

So the quenched polymer distribution Qu,v is the distribution of the backward Markov 
chain with initial state v, transition and absorption at u. The distributions Qu,v 
are the noisy counterparts of geodesies. The nesting property of geodesies manifests itself 
through conditioning. Let u < z < w < v in Z^. Let Az^w be the set of paths in IIu^^ that 
go through the points z and w. Given y_ G u,, let By be the set of paths in n„ „ that 
traverse the path y, (that is, contain y. as a subpath). Then 

(2.9) QuA^y. I = QzAv.}- 

Define the random geodesic spanning tree Tu rooted at u by choosing, for each x G 
(n + Z^) \ {u}, a parent 

ei with probability "^"(x,x — ei) 
62 with probability '^"(x,x — 62). 



(2.10) 7(x) 



Now that we have the positive temperature counterparts of geodesies, we can find the 
positive temperature counterpart of the competition interface, by reference to the tree To 
rooted at 0. Let 7o,er- be the subtree rooted at e^, so that To = {0} U To,ei U To,e2 as a 
disjoint union. The lemma below shows that there is a well-defined path X that separates 
the trees To,ei and To, 62 > and evolves as a Markov chain in the environment defined by the 
partition functions. In other words, this random walk in a random environment (RWRE) 
is the positive temperature analogue of the competition interface. The picture for X_ is the 
same as for ip in Figure 1. 
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Lemma 2.2. (a) Given the choices made in (2.10), there is a unique lattice path {Xn)nez+ 
with these properties: Xq = 0, X„ — Xn-i G {61,62}, and for each n and r € {1,2}, 
{X,^ + ker : A; G N} C To,er- 

(b) Xn is a Markov chain with transition matrix 

(2.11) 7r,,,+e. = 1 "'T":",! , X G Z^, r e {1, 2}. 

Proof, (a) To prove the existence of the path, start with Xq = 0, and iterate the following 
move: if j{Xn + ei + 62) = X„ + 6^ set Xn+i = Xn + e^-r- 

(b) Given the path (^fc)fc=o with Xn = x, we choose Xn+i = x+ei if 7(x+ei+e2) = x+e2 
which happens with probability 

Similarly for Xn+i = x + 62 with the complementary probability. □ 

A genuine RWRE transition probability should satisfy Trx,y{uj) = T:Q^y-x{TxOj) for shift 
mappings (Ta;)^g^2^ acting on the environments uj. We augment the space of weights to 
achieve this. We need to be precise about the sets of sites on which various classes of 
weights are defined. 

Definition 2.3. A collection of positive real weights (^,fy,C)0 = {Cx, f/x-e2) Cx-en Ca;-ei-e2 ■ 
X € N^} satisfies north-east (NE) induction if these equations hold for each x € N^: 

(2.12) Vx = Cx , Cx = i ^""-'^ 



I > 5 SX <,X f. ) 

^x— 62 ~r Sx—ei 'nx—e2 ' ei 

(2.13) and ei— 62 = ^x— £2 + Cx—ei- 



North-east induction simply keeps track of ratios of partition functions. Take as given 
the subcollection of weights {CijiVifiiCoj '■ hj ^ I^} on and construct the polymer 
partition functions 



(2.14) 

Then define 
(2.15) 



0,1) 



Ml 

^ l[V-' with Vij = < 

X. Gilo, I, i=l 



(i,j) €N2 
i G N,j = 
i = oj G N. 



^O.x— ei 



and 



-'0,x—e2 



for X e N^. 



^0,x ^0,x 

Now the subsystem (^,r/, C) satisfies (2.12), as can be verified by induction. To get the full 
system (C, C, just define 4 = Vx+e^ + Cx+ea for x £Z\. 

Note that (2.15) is valid also on the boundaries, by the definition (2.14) of Z^^kcr for 
/c € N. The reason for the distinct notation {f?i,0)Co,i} for boundary weights in (2.14) is 
that these are also ratios of partition functions, just as r]x and C,x in (2.15). We shall find 
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that in the interesting log-gamma models the boundary weights {rji^Q, Co,j} are different from 
the bulk weights {dj}- The role of the ^ weights is not clear yet, but they will become 
central in the log-gamma context. 
Define the space of environments 

O / fC ^ ^N2 + (NxZ+)+(Z+xN)+Z^ 

^2 = {u} = V, C, e • 

(^)^)C)0 satisfies NE induction}. 

Translations act naturally: for z G Z^, T^u = (^2+n2, f/2+NxZ+, Cz+z+xN, where we 

introduced notation i^2-|_n2 = {S,z+x} xeN'^ 7 similarly for the other configurations. 

Definition 2.4. The polymer random walk in random environment is a RWRE with envi- 
ronment space One and transition probability 

(2.17) ■n-x,x+ei{^) = ^^^7^7 and ir^^^+eii^) ~ '^'"'^''^ 



u ■■ j,,j,-ro^ \— / L. 

This definition is the same as (2.11) with partition functions (2.14). The quenched path 
probabilities P'^ of this RWRE started at are defined by 

n 

(2.18) P'^{Xo = 0,Xi=Xi,...,Xn=Xn) = '[[7rx,_„xM withxo = 0. 

k=l 

Distributions of are again related to polymer distributions and partition functions. 
Define 

(2.19) Zo,„ = 5^ n ^ ^ 

a;.eno,u i=0 

In contrast with (2.14), this time the weight at the origin is included but the weight at v 
excluded. From (2.13) and (2.18) we derive two formulas. First, distribution of X„ is a 
ratio of partition functions: 

(2.20) p^{Xn = x) = ^ for X G such that = n. 

Then, if the walk is conditioned to go through a point, the distribution of the path segment 
is the polymer probability in ^ weights: for x, = (xfc)^^Q € ^o,xn, 

^ n—l 

(2.21) P^iXo = 0,Xi=Xi,...,Xn = Xn\Xn=Xn) = J J] = ^O.xJ^;.}. 

3. The log-gamma polymer 

This section gives a quick definition of the log-gamma polymer and its Burke property. 
Let < A < p < oo. 

Definition 3.1. A collection (C,r/,C,I) = {Cx, ■??x-e2, Cx-ei, Ix-ei-ea : x € N^} of positive 
real random variables is a gamma system of weights with parameters (A, p) if the following 
three properties hold. 
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(a) NE induction (Definition 2.3) holds: for each a; G N^, almost surely, 



(3.1) 



Vx—e2 



Cx — Cx 



Cx—ei 



Vx—e2 ' >,x—ei Vx—e2 ' ^x—ei 

and ^x—ei—e2 ~ 'nx—e2 ~l~ Ca;— ei- 

(b) The marginal distributions of the variables are 

(3.2) rjx ~ Gamma(A), ~ Gamma(p — A), and ~ Gamma(p). 

(c) The variables {CijiVifiXoj ■ hJ ^ are mutually independent. 

A triple rj, is a gamma system with parameters (A, p) if conditions (a)-(c) are satisfied 
without the conditions on 

Note that the ^-weights are defined only in the bulk N^, while the ^-weights are defined 
also on the boundaries and at the origin. Variables rj^ and C,x can be thought of as weights 
on the edges, while and Cx are weights on the vertices. Edge weights will also be denoted 
by 

(3.3) Tx—ei,x — f]x and Tx—e2,x — Cx- 

A natural way to think about equations (3.1) for a fixed x is as a mapping of triples 
(Figure 2): 

(3.4) (Cx) ??a;— 62 ) Cx-ei ) (??x , Cx ) Cx-ei-e2 ) • 



5X— ei 




Figure 2. Mapping (3.4) that involves variables on a single lattice square. The 
picture illustrates how southwest corners are flipped into northeast corners in an 
inductive proof of the Burke property. 

This mapping preserves the independent gamma distributions: if (Cx, ^x-e2) Cx-ei) are 
independent with marginals (3.2), then the same is true for (r/^;, Cxj Cx-ei-e2)! as can be 
checked for example via Laplace transforms. Consequently a gamma system (C) C; C) can 
be constructed by repeated application of equations (3.1) to independent gamma variables 
given in (c). 

An equivalent way to define the gamma system is to first construct the following polymer 
partition functions from the weights given in (c): for < u < v in Z^, 



(3.5) 



■'u,v 



x.€Tlu,v i=l 



lCo,i, i = 0,jGN. 
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Then define 

(3.6) Tx-er,x = — 1^—^ for X € and r G {1, 2} such that x — e,. S Z^. 

The weights rjx and are then defined via (3.3). Now we have a gamma system (CvX): 
which can be augmented to a gamma system (^, r/, ^) since ^ is a function of {r], (")• 

Mapping (3.4) furnishes the induction step in the proof of the Burke property of the 
log-gamma polymer [27, Thm. 3.3]: for any down-right path on Z^, the r- variables on the 
path, the ^ variables strictly to the northeast of the path, and the ^ variables strictly to the 
southwest of the path are all mutually independent with marginal distributions (3.2). The 
induction proof begins with the path that consists of the ei- and e2-axes. Southwest corners 
of the path can be flipped into northeast corners by an application of (3.4), as illustrated 
in Figure 2. 

As an application of the Burke property consider the down-right path consisting of the 
north and east boundaries of the rectangle {0, . . . , m} x {0, . . . , n}. Then the Burke property 
gives us this statement: 

variables {rji^n, Cm,j, li-ij-i ■ I < i < m,l < j < n} 
are mutually independent with marginals (3.2). 

Remark 3.2. Let us revisit the zero temperature limit (Remark 2.1). The log-gamma 
polymer does not have an explicit /3 parameter but p represents temperature. Replace p 
by ep in the definitions above, so that ^a; ~ Gamma(e/9). Then as e \ 0, —elogS^x =^ ^^x, 
a rate p exponential weight. For u < v in N, 

\v-u\i 

e log = e log ^ expj -e"^ ^ elog^^^j 

X.€Tlu,v «=1 

1)— n| 1 



ujx = Gu,v as e \ 0. 



^ max 

X,^^u,v 

1=1 

In other words, we have convergence in distribution to last-passage percolation with expo- 
nential weights. 

An important function of the polymer path is the exit point or exit time toxit of the path 
from the boundary: t^xit = tei V , 

(3.8) te^ = max{A; > : = (i, 0) for < i < k] 
and 

(3.9) iea = max{£ > : = (0, j) for < j < £}. 

Note that for each path te^ = 0. Partition functions (3.5) based at can be equivalently 
written as 

*oxit I I'll 

(3.10) zo,.= \Vl<-^.^))\ n ^^.0' ^^^+- 

x.&Iq„ i=\ i=*exit + l 
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In a {X, p) gamma system we have the means E(log7yj^o) = ^o(A) and 
(3.11) E(log Zo,(„,„)) = -m^o(A) - n^o{p - A). 

The second one comes from 



n 

(3-12) logZo^(„^„) = -^logrji^o -^logCm,j, 

1=1 j=i 

a sum of two correlated sums of i.i.d. random variables. Above = T'/F is the digamma 
function. It is strictly increasing on (0,oo), with \I'o(0+) = — oo and ^'o(oo) = oo. Its 
derivative is the trigamma function ^'i = ^'q that is convex, strictly decreasing, with 
^'i(0+) = oo and ^i(cx)) = 0. 

The asymptotic directions (or velocities) of admissible paths in lie in the simplex 
U = {u = {u,l — u) : u G [0, 1]}. Fundamental to the behavior of the log-gamma polymer 
is a 1-1 correspondence between velocities u £U and parameters A G [0,p], for a fixed p. 
The characteristic direction for (A, p) is 

^i(p-A) ^i(A) 

,^i(A) + ^'i(p- A)' ^'i(A) + ^i(p- A)- 

Conversely, for u = (u, 1 — u), the unique parameter 6{u) = 9{u) € [0, p] for which u is the 
characteristic direction is defined by 6{0) = 0, 6{1) = p, and 

(3.14) -u^i{e{u)) + {l-u)^i{p-e{u)) = for nG (0,1). 

Function 6{u) is a strictly increasing bijective mapping between u € [0, 1] and 6 G [0,p]. 

The function 9{u) will appear throughout the paper. Let us point out that if {ni,n) = 
cu then the right-hand side of (3.11) is minimized by A = 9{u). As we shall see, this 
identifies the limiting free energy for the log-gamma polymer with i.i.d. Gamma(/9) weights. 
Notationally, A, a, v denote generic parameters in [0, /o], while is reserved for the function 
defined above. 

4. Limits of ratios of point-to-point partition functions 

Fix < /5 < oo. Let i.i.d. Gamma(p) weights w = {wx ■ x € Z^} be given. Define 
partition functions 

\v—u\i — l 

(4,1) Z„= n '><'"< "in 



'^Xi 1 U ^ U ^ (7 Hi Z^ 

I'.eriu,!, i=o 



Note that the weight at u is included and v excluded, in contrast with definitions (2.5) and 
(3.5). This is for convenience, to have clean limit statements below. 

Suppose a lattice point (m, n) G tends to infinity in the first quadrant so that it has 
an asymptotic direction in the interior of the quadrant. Let A G (0, p) be the unique value 
such that this assumption holds: 

r 4 1 "T- ^'l(P — A) 

(4.2) mAn-)-oo and > — . 

n ^'i(A) 

When (4.2) holds we say that {m,n) — )• oo in the characteristic direction of {X,p)- 
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The central theorem of this paper constructs gamma systems out of i.i.d. weights by 
taking hmits of ratios of point-to-point partition functions. 

Theorem 4.1. On the probability space of the i.i.d. Gamma(p) weights w = {wx ■ x G 
Z^} there exist random variables {C^, r/^.gj, C^-ei • ^ (0, p), x € N^} with the following 
properties. 

(i) For each X G (0,p), {^■'^ , r]'^ , ^ w) is a gamma system with parameters {X,p). Fur- 
thermore, if on the same probability space there are additional random variables rj, () = 
{■^xj f?x-e2) Cx-ei '■ X G N^} such that iCrjX^w) is a gamma system with parameters {i^,p), 
then (C,^,C) = (e^^/^C) a.s. 

(ii) If a sequence {m,n) — >■ oo m the characteristic direction of{X,p), as defined in (4.2), 
then these almost sure limits hold: 

(m,n)-^oo ^x--ei,(m,n) {m,n)-j.oo ^jy_e2,(m,n) 

for a// X E N X Z+ and ?/ € Z+ x N. 

(iii) The weights are continuous in X, and the edge weights are monotone in X: for each 
X for which the weights are defined, almost surely, 

(4.4) vt<V^' and Ct>C^' forXi<X2 
and 

(4.5) rj'x ^ ^ C, ^ C asX^u. 



The rest of this section proves Theorem 4.1. The reader not interested in the (rather 
technical) proof can proceed to the next section where these limits are applied to solve a 
variational problem for the limiting free energy. 

The proof relies on the following lemma for gamma systems. Let (^, 77, ^) be an (a, p)- 
system according to Definition 3.1. Using the ^ weights, define partition functions 

\v--u\\ — l 

(4.6) Z,,, = n (^'-0^'' 0<n<t;inZ^, 
and for 2; € N x Z+ and y € x N edge ratio weights 

(4.7) Ix,{m,n) = ' ' ^-Ild Jy^[m,n) = ^ ^ ■ 

•^x— ei ,(m,n) ^y—e2,(m,n) 

Lemma 4.2. Lei 0<A<a<A<p. Consider two sequences {mi,ni) — )• 00 and {fhj,nj) — >■ 
00 in such that 

(4.8) m,^5;i(p_A) m^^ Mp-~^) _ 

rii ^'i(A) nj ^'i(A) 

Then /or x G N x Z+ and y G Z+ x N 

(4-9) lim 4,(m„n,) < f?x < lim 4,(mj,nj) a-s- 
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and 

(4.10) lim Jj/,(^,,nj) < Cy < lim Jy,im„n,) a.s. 



Proof of Lemma 4- ^- For notational simplicity we drop the i,j indices from {m,n) and 
(m, n). We relate ratios (4.7) to ratios of partition functions with boundaries. Let ^^^^-^ 
denote a partition function that uses r] and C weights on the north and east boundaries of 
the rectangle {k, . . . , m} x {£,..., n} and ^ weights in the bulk: 

m ^ " 1 

^(k,n),(m,n) = H ~ ' ^{m,e),(m,n) = Jl 7 

s=k+l <=£+l 

and for < k < m and < ^ < n 

m— 1 

7NE 

^(fc,^),(m,n) ~ / > c J.X n 

(4.11) 



n— 1 ^ n ^ 

+ ^^[k/)Xra~-i,jn — : n 7 — • 



In the last formula Z'^^ n) ^® decomposed according to the entry points (i, n) and {m,j) 
of the paths on the north and east boundaries. If the entry is at {i,n), the first boundary 
variable encountered is ?/i+i,n associated to the edge {{i, n), (i + 1, n)}. The last bulk weight 
S,i,n~i has to be inserted explicitly into the formula because ^(fc/),(i,n) does not include the 
weight at (i,n), by its definition (4.6). 

The corresponding ratio weights on edges are 

. . (fc,£),(m,n) , (fc/),(m,n) 

(4.12) ^(fc,Q,(m,n) = and J{k,e),im,n) = 1^ • 

^{k-l/),{m,n) ^(fc/-l),(m,n) 

Due to the reversibility of the shift-invariant setting, these ratio weights are the same as 
the original ratio weights, and thereby do not depend on (m,n). This is the content of the 
next lemma. 

Lemma 4.3. For < k < m and < i < n such that the weights below are defined, 

(4-13) "Hk/ = I(k,e),{m,n) 0.nd Ck/ = J(k,e),{m,n)- 



Proof of Lemma 4-3- When i = n in the r/-identity or A; = m in the (^-identity the claims 
follow from the definitions. Here is the induction step for r/^. £, assuming the identities have 
been verified for the edges {{k — 1,£ + 1), {k,i + 1)} and {{k,i), {k,£ + 1)}, closest to the 
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north and east of the edge {{k — {k,£)}: 



7NE 



^ _ ^fc-l/^(fc/),(m,re) 

^{k,l),{m,n) r^NE _|_ r^NE 

{k,e),{m,n) (A;-l/+l), 
"^(A.-l,^+l),(m,n) _ ^(fc/+l),(m,n) ^ 

' r^NE ^NE 

^(fc,^+l),(m,n) ^{fc,^),{m,n) 

= 4-i,£(l + ] = (%,^ + Cfc-i,^+i)( 1 + ^tlM±l\ = r]f,i. 

V / V Vk,i J 

The third equahty is the induction step. The fourth equahty uses (3.1) twice. □ 

We need one more variant of ratio weights, namely the types where the last step of the 
path is restricted to either ei or 62- Relative to any fixed rectangle {k, . . . , m} x {£,... , n}, 
define the distances of the entrance points of the polymer path x, G 'n^{k,e),{m,n) on the north 
and east boundaries to the corner {m,n): 

(4.14) t*^ = max{r > : Xm-k+n-e-i = (m - i, n) for < i < r} 
and 

(4.15) = max{r > : Xm-k+n-l-j = (m, n - j) for < j < r}. 

For a subset A of paths, write Z{A) for the partition function of paths restricted to A (in 
other words, for the unnormalized polymer measure). Then define, for r € {1,2}, 

^ ^(M),(m,«)(^er > Q) 
{k,e),im,n) y 0) 

_ ^(fc/),(m,n)lV > 

We are ready to prove Lemma 4.2. We go through the proof of (4.9), the case for J 
being the same. Applying Lemma 9.1 from the appendix (to a reversed rectangle) gives 

(4.17) %,£ + - %,^) < - ^fc/ + (■^(M),(m+l,n+l) - • 

Taking (4.13) into consideration, the task is 

< < \^ ~ -^(fc,^),(m+l,H+l)}- 

(m,n)— ^00 

We do the first limit for 62. The second is similar. Introduce the parameter 

^- i,.(,-A) + i.i(A) ^°° 4r^-(*,(p-A),*,(A)). 

The first inequality of (4.18) follows from showing that Ve > 3a > such that 

(4-19) nilk,eUm+l,n+l) > Iik,eUm+l,n+l) + s} < 26"'^^. 
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Introduce the quenched path measure Q^^ ^-^ „) that corresponds to the partition func- 
tion in (4.11). 

62 _ ^(fe/),(m+l,n+l)(^e2 > 0) 

(fc,^),(m+l,n+l) ~ Ttne 77* ^Tvi 7"^ 

^(fc-l,£),(m+l,n+l)'^''e2 ^ ^ (k-l/),{m+l,n+l) 



< 



(fc/),(m+l,n+l) 



-l,^),(m+l,n+l) 

For small enough e^v the probability in (4.19) is bounded above by the sum 

(4.20) > ^} +F{Qlk-i,iUm+i,n+i)K > 0) > en}. 

Note that the event in the Q""^ -probability was replaced by its complement. A sequence 
< etv \ will be chosen below. 

By (4.13) I(k,e),{m+i,n+i) has Gamma(a) distribution and so the first probability in (4.20) 
is bounded by e~"^l'^'^ . 

We show that the Q'^^ -probability in (4.20) is actually a large deviation by replacing 
{m,n) with a direction that is characteristic for (a,p). The next lemma contains the idea 
for replacing (m,n). 

Lemma 4.4. Let {fh,n) satisfy rh> m and i < n < n. Then 

(4-21) Qlk/Um,n)it*ei > 0) = > rfi - m). 

Proof. A path in Il(k,e),{m,n) that satisfies t*^ > must use one of the edges {{i,n — l), {i, n)}, 
k < i < m—1. Otherwise it hits the east boundary first and t*^ = 0. Decomposing according 
to this choice of edge and using definition (4.11), 

■m— 1 ^ 



i=k ^'{k,i),{m,n) 

By Lemma 4.3 the last ratio does not depend on {m,n), and {m,n) can be replaced by 
(fh,n). This moves the northeast corner in definition (4.11) to {fh,n), as well as the 
reference point of t*^ in (4.14). Since the sum still runs up to m — 1, it now represents 
paths in 'n^(k.e),{m,n) that hit the north boundary to the left of (m, n). This proves Lemma 
4.4. ' ' ' □ 

Take (fh,n) = ( [A^^i(/9 — a)J + k — 1, [A^^'i(a)J + £), essentially the characteristic 
direction for {a,p). Since A < a and ^'i is strictly decreasing, there exists 7 > such that 
for large enough N, fh > m + l + Nj and n < n — N"f. Put sn = e~^^'^^ for a small enough 
5i > 0. Then for large enough N ^ 

IP{Q(fc-i,.),(m+i,n+i)(*:. > 0) > em] < P{Q(^i,,),(^,^)(4 > Nj) > e-^^^^} 

The last inequality came from Lemma 9.2 in the appendix where we can take kn = 1 and 
(5 < 7. Both probabilities in (4.20) have been shown to decay exponentially in N, and 
consequently (4.19) holds. This completes the proof of Lemma 4.2. □ 
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Turning to the proof of Theorem 4.1, we begin by showing the a.s. convergence in (4.3) 
for a fixed sequence and fixed A. Define ratio variables by 

/ , „qN ^x,{m,n) A /- ^yA^>^) 

•^x— ei,(m,n) ^y—e2,{m,n) 

for X G N X Z+ and y G Z+ x N. 

Proposition 4.5. Fix < A < p and fix a sequence {m,n) oo as in (4.2). Then for all 
a; G N X Z+ and y G x N the almost sure limits 

(4.24) r]^ = hm Vx,(m,n) and Cy = , 1™ ??j;,(m,n) 

(m,n.)— >oo (m,n)— >oo 

exist and have distributions rjx ~ Gamma{X) and C,y ~ Gamma{p — A) . 

Proof. We treat the case of the r] variables, the case for ^ being identical. For a while, 
until otherwise indicated, we are considering a fixed sequence of lattice points that satisfies 
(m, n) ^ oo as in (4.2). To avoid extra notation we refrain from indexing the lattice points, 
as in (mfc,nfc). Later we can improve the result so that the limit only depends on A and 
not on the particular sequence (m,n) oo. 

We show that for < s < oo the distribution functions 

(4.25) G*(s)=P{ lim VxJm,n) < s} and G*(s)=P{ lim r]x(m,n)<s} 

{m,n)-i-oo (m,n)-j.oo 

satisfy G*{s) = G*(s) = Fx{s) where 

Fx{s)=r{X)-^ f t^~^e-Ut 
Jo 

is the c.d.f. of the Gamma(A) distribution. Since lim y^. jm r,) < lim?]^. this suffices for 
the conclusion. Working with the distributions allows us to use any particular construction 
of the processes. 

Let {Uij} be i.i.d. Uniform(0, 1) random variables. For i,j G N and a G (0,/>) define 

(4.26) f]t^^ = F-\U,,o) and CS, = F;\{Uo,j). 

This gives us coupled weights fjf^ ~ Gamma(a) on the south boundary and Co j ~ Gamma(yO— 
a) on the west boundary of the positive quadrant. For the bulk weights take an i.i.d. col- 
lection {crx}x^n^ of Gamma(/9) weights independent of {Uij}. 

As mentioned after Definition 3.1, the mutually independent initial weights {o"jj-, f/fg, Cq j : 
i,j G N} can be extended to the full gamma (a, p) system (cr, C", o"^"^)- The construction 
preserves monotonicity of the edge weights, so that 

(4.27) f?-. < f^lj and Q > CI, for a<v. 

Bracketed superscript [a] reminds us that even though the variables {o'Yj}ij>o are i.i.d. 
Gamma(/9) for each a G (0,yo), they were computed from a-boundary conditions. Define 
partition functions 

•u— 'U| 1 — 1 

(4.28) E n Q<u<vinZ\ 
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and edge ratio weights 

(4.29) P] ,= "'('"'"^ and v^^^J^^ 

x—ei,{m,n) y—e2,{m,n) 

For each a G (0, p), we have equahty in distribution of processes: 

(4.30) {^{i+l,j),{m,n)^ ■^{ij+l),(m,n)^^lj} = {'n{i+l,j),{m,n) X{i,j+l),{m,n) , Wij} . 

These processes are indexed by {(i,j) G Z^, {m,n) > (i + l,j + 1)}. The equahty in 
distribution comes from identical constructions apphed to i.i.d. Gamma(p) weights: on the 

left to at"!, on the right to w. Now in (4.25) we can use any process 

Applying Lemma 4.2 to two gamma systems (cr, r/"^ , (^"^ , (t'"^!) and (cr, r/"^ , C"^, (t'^^I) 
gives, for any < ai < A < 02 < /O, 

(4-31) lim /[-l > and lim < r^- a.s. 

(m,n)— >oo (,m,?ij— >cx3 

By the equality in distribution (4.30), 

G,(s)=P{ hm l];^%^^^^^<s]<V{fl'^^,<s} = F^M\Fx{s) as ai / A, 

(m,n)—)- 00 

and 

G*(s) = P{ lEI /["^l . . < 4 > F,,(s) Fa(5) as a2 \ A. 

(m,n)— >oo ^ ' ' ' 

This gives Fx{s) < G*{s) < G*(s) < Fx{s) and completes the proof of Proposition 4.5. □ 

Proposition 4.5 gave the a.s. convergence of ratios along a fixed sequence and for a given 
A € (0, /o). Next we construct a system of weights (^,r/, ^,?i;) from the limits (4.24) by 
defining 

(4.32) ^. = Vx + Cx, forxGN^. 

Proposition 4.6. The collection (^, tj, w) is a gamma system of weights with parameters 
(A,p), that is, it satisfies Definition 3.1. 

Proof. Equations (3.1) follow from the limits (4.24) and 

^x+eiXfn.n) ~l~ ^x+e2,{m,n) 

= ^ ■ 

^x,{m,n) 

By the equality in distribution in (4.30), it also follows that the limits in (4.31) exist: 
Let < ai < A < Q!2 < p. Utilizing (4.30), (4.31), and (4.33), 

irj,w) ^ (/M,d'"^i) > (rs<^'°^i) — (^\<?t") 

(4.34) 

and (r/,u;) = (/["2l,a["2l) < (j?°2,(7["2]) _^ {f]^,a^^^). 

02\A 
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The inequalities and the convergence are a.s. and coordinatewise. The convergence follows 
from the continuity of the definitions (4.26) in a and the continuity in the equations (3.1) 
that inductively define the (f/"^, C", d"'"]) weights. The consequence is that 

(4.35) (r/,u;)^(f/\dW). 

Equations C,x = Wx-e2 ~ ^x-e2+ei and = Vx + Cx map (ry, w) to the full system 
{S^jTiXtw). The same mapping applied to the right-hand side of (4.35) recreates the sys- 
tem (cj, f]'^, C'^, (jl'^l), which we know to be a (A, p) gamma system by its construction below 
(4.26). □ 

Proof of Theorem 4-1- Fix a countable dense subset D of {0,p) and VA € -D a sequence 
{m,n) oo that satisfies (4.2). By Propositions 4.5 and 4.6, we can use limits (4.3) along 
these particular sequences to define, almost surely, (A, p) gamma systems (^^, r/^, (^'^, tw) for 
A G L*. Monotonicity (4.4) is satisfied a.s. for Ai, A2 G -D by Lemma 4.2. (The point is that 
the Z partition functions in (4.6) are the same for all systems r]'^ , , w) .) 

Monotonicity and known gamma distributions give also the limits in (4.5) when A — )• 1/ 
in D. For example, suppose A u in D. Then limx/'ur]^ < r]^, but both are Gamma(z^) 
distributed and hence coincide a.s. The limit — )• comes from the limits of rj and C and 

^x Vx ~^ • 

Extend the weights to all AG (0, p) by defining 

(4.36) rj^ = inf{r?,^ : u G D n {X, p)} = snp{r]^ :aGDn{0,X)} 

with the obvious counterpart for (^^ and then — Vx ~^ Cx- The inf and the sup in (4.36) 
must agree a.s. because (i) the sup is not above the inf on account of the monotonicity for 
X £ D, and (ii) they are both Gamma(A) distributed. By the same reasoning, for A G -D 
definition (4.36) gives a.s. back the same value Vx originally constructed. 

To check that the new system {C'^ , v^ ^ : is a (A,p) gamma system, fix a sequence 
D B Ui X, and observe that equations (3.1) are preserved by limits and the correct 
distributions come also through the limit. Extending properties (iii) utilizes monotonicity 
again. Limits (4.3) of ratios for a value X ^ D come from Lemma 4.2. 

The final point is the uniqueness in part (i). Lemma 4.2 and the limits (4.3) imply that 
Vx^ < ^fe < Vx^ for sll ai < u < a2- □ 



5. BUSEMANN FUNCTIONS AND A VARIATIONAL CHARACTERIZATION OF THE FREE 

ENERGY 

In this section we turn the limits of ratios of point-to-point partition functions into 
Busemann functions, and use these to solve a variational formula for the limiting free energy. 
The parts from this section needed for the sequel are definition (5.12) of the velocity u(/i) 
associated to a tilt h, and the large deviation bound (5.20). The latter is needed for the 
proofs in Section 6. 

We consider briefly general i.i.d. weights w = {wx)x^i?^ on a probability space (0,(3,P) 
assumed to satisfy 

(5.1) 3e>0: E( [ log u;oP+^) < 00. 
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Later we specialize back to wq ~ Gamma(p) . It is convenient to use exponential Boltzmann- 
Gibbs factors. Let p{ei) = p{e2) = 1/2 be the kernel of the background random walk Xn 
with expectation E and initial point Xq = 0. Define the potential g{w) = — log wq + log 2. 
In this notation the point-to-point partition function (4.1) is 

(5.2) Zo,, = E[e^"^=o9iTx,^),Xn = v], n = \v\i. 
Introduce a tilted point-to-line partition function 

(5.3) Zo^^(^) = E[e^k=o 9iTx^^)+h■Xt,^^ ^ ^ (/ii, /12) G and iV G N. 

The set of limit velocities for admissible walks in \s Li = {(n, 1 — u) : < n < 1}, 
with relative interior milA = {(n, 1 — ti) : < n < 1}. For each u = (ti, 1 — n) € Z//, let 
Xn(u) = (LnuJ,n — [ntij). Define limiting point-to-point free energies 

(5.4) Ap2p(u) = lim n"MogZo^„(u), u^U, 
and tilted point-to-line free energies 

(5.5) Ap2,(/i) = lim N-^ log Z^.^., h = (/ii, /la) G M^. 

Under assumption (5.1) these limits exist P-a.s., Ap2p is continuous and concave in u, and 
Ap2£ is continuous and convex in h [24]. 

We recall two variational formulas, valid for i.i.d. weights under assumption (5.1). First, 
a convex duality between the free energies ([24, Rmk. 4.2], also proved below in (5.19)): 

(5.6) Ap2p(u) = inf {kp2e.{h) -u-h}. 

Let /C denote the class of correctors F : Qx {ei, 62} — )• K that satisfy F S L^, E,F{w, z) = 
for z £ {61,62}, and a cocycle property: F{w,ei) + F{Tej^'w,e2) = F{w,e2) + F{Te2W,ei) 
P-a.s. Then we have the variational formula [25, Thm. 2.3] 

(5.7) Ap2£(/i) = inf P-ess sup log V p(^z)e3M+''-'+^^'"''\ 

F&K. w f — ' 

We solve (5.6) and (5.7) for the log-gamma model. The next corollary turns the limits of 
Theorem 4.1 into Busemann functions, and states the properties needed for the development 
that follows. Recall the function 6{u) € [0,p] of (3.14), the unique parameter such that u 
is the characteristic direction for {6{u),p). 

Corollary 5.1 (Corollary of Theorem 4.1). Assume {wx} are i.i.d. Gamma{p). 

(a) For each velocity u € intU and for each x,v £ the F -almost sure limit 

(5.8) B'^{w,x)= lim (logZo,^,„(u)+^ - logZa.^„(u)+„) 

exists and is independent of v. 

(b) The sequences {B^ {Ti^j^w , ei) : i € Z+j and {B^(Tji.2W,e2) ■ j G Z+} are i.i.d. with 
g-B"{«;,ei) _ Gamma{0{u)) and e-^"^'"'"^) ^ Gammaip - 0{u)). 
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We call i?" a Busemann function, by analogy with the Busemann functions of last- 
passage percolation which would be limits of differences Go^i„(u)+v — Gx^xn{u.)+v Of course 
we are merely re-expressing limits (4.3) in the form 

„-B"(T,u>,ei) _ 1- -^x+ei,x„(u)+i» _ e(u) 
— nm — '/x+ei- 

The admission of the perturbation v in (5.8) gives the cocycle property: 

(5.9) B'^iw, x) + B^'iT^w, y) = B'^iw, x + y). 
As a function of u € int^, define the tilt vector 

2 

(5.10) /i(u) = (/ii(u),/i2(u)) =-J]E[B"(u;,ei)]ei = (^'o(e(u)),^'o(p-^(u))). 

i=l 

Note that /i(u) is not well-defined for u on the axes, ^(u) converges to (to p) as u 
approaches the y-axis (x-axis). Then one of the coordinates of /i(u) approaches — oo. The 
function 

(5.11) u = (n, 1 - u) ^ /ii(u) - /i2(u) = ^'o(0(u)) - ^'o(p - 0(u)) 

is a continuous, strictly increasing function from u G (0,1) onto ( — 00,00). An inverse 
function to (5.10), 9 /i i-t- u(/i) G milA, is given by 

u = u(/i) uniquely characterized by the equation 

^^'^'^^ hi-h2 = MOH) - -^oip - 0{u)). 

Note that u(/i) is constant when h ranges along a 45 degree diagonal. If /i = there is no 
tilt, u(0) = (1/2, 1/2), and 6l(u(0)) = p/2. 
From these ingredients we solve (5.6). 

Theorem 5.2. Let u = (u, 1 — n) € intiY. TiZi /i(u) kills the point-to-line free energy: 
Ap2£(/i(u)) = Vu € int^. Furthermore, h{u) minimizes in (5.6) and so 

(5.13) Ap2p(u) = -u • h{u) = -u^oieiu)) - (1 - n)^o(p - ^(u)). 



Define the corrector 

(5.14) F"(u;,z) = -5"(ii;,z) -/i(u) -z, zG {61,63}. 
Theorem 5.3. Given h = {hi, /12) G M^, t/ie equation 

(5.15) /ii(u) - /i2(u) = /li - /l2 

determines a unique u G intZ//. T/ien ^ K, is a minimizer in (5.7). T/ie right-hand side 
of (5.7) is constant in w so the essential supremum can be dropped: ¥-a.s., 

A,2e{h) = log V p(z)e9(-)+'^-+^"(-'^) = -h^iu) + h^ 

(5.16) zen 

= -*o(p-^(u)) + /l2- 
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Remark 5.4. Theorem 5.2 is the third proof of the exphcit value of Ap2p(u). This result was 
first derived in [27] together with fluctuation bounds. The simplest proof is in [15] where 
the minimization of the limit of the right-hand side of (3.11) is done with convex analysis. 
The value (5.16) of the tilted point-to-line free energy has not been computed before. 

Remark 5.5 (Large deviations). Let us observe how the duality between tilt h and velocity 
u in (5.6) is a standard large deviation duality. The tilted quenched path measure is 

(5.17) <(iV)K} = -^eS£o^^(^^.-)+'^-^-P{x.}. 

0,(Af) 

The quenched large deviation rate function for the velocity is 

A(v) = - lim M logQ^^J \N-^Xn - v| < 5} (P-a.s.) 

(5.18) N^oo 

= Kp2i{h) - h-v - Ap2p(v). 

The last equality uses the continuity of Ap2p and Lemma 2.9 in [24]. The limiting logarithmic 
moment generating function is 

N^oo 

By Varadhan's theorem these are convex duals of each other: 



AQ,hia) = lim N'^ log E<W [e"-^^] = Ap2,{h + a) - Ap2,{h) P-a.s. 



(5.19) A(v) = sup {a • V - AQ,h{a)} 

agR2 

which is the same as (5.6). For the next section we need the minimizer of Ih- By (3.14), 
(5.13) and calculus, Ih is uniquely minimized by u(/i) defined by (5.12). Consequently the 
walk converges exponentially fast: for 6 > 0, 

(5.20) IhR N-^logQo(^.{\N-^XN -u{h)\>6} <0 F-a.s. 

The function Ap2p(u) extends naturally to all of by homogeneity: Ap2p(cu) = cAp2p(u). 
Part of the duality setting is that the mean of the Busemann function gives the gradient 
VuAp2p(u) = -/i(u). 

The remainder of this section proves the theorems. 

Proof of Theorem 5.2. That is a corrector is clear by (5.9). Let 

m— 1 

/"(w, x) = ^ F'^{T^^w, Xi+i - Xi) = -B'^iw, x) - h{u) ■ x 

i=0 

be the path integral of F. The admissible path {xjj^o above satisfies xq = and Xm = x, 
and the cocycle property implies that /" depends on the path only through the endpoint 
X. Corollary 5.1(b) verifies exactly the sufficient condition (9.11) for (9.10). From Theorem 
9.3 in the appendix, 

\f''iw,x)\ 

(5.21) max '-^ ^ ' ^' a.s. 

x€Z'l:\x\i=n n 
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This ergodic property slips /"(w, Xn) into the exponent in the free energy hmit, and shows 
that tilt h{u) kills the point-to-line free energy: 

71— >00 

(5.22) n^oo ^ ^ J 

= lim n-hogE[e^"^=oi9(Tx,w)+h{uyz,+,+F-{T^^w,z,+,)h 

n— !>oo 

= 0. 

The last equality comes from 



(5-23) _ ^^p(z)ef(-)Z 



-2'o,ai„(u) -^0,£„(u) 



Fix xi ^lA. Since \Xn\i = n, the expression on the right-hand side of (5.6) satisfies 

K2i{h) -u-h = Ap2iihi - /i2, 0) - u • {hi - /i2, 0) 

and so, as a function of h, is constant along 45 degree diagonals. So the minimization needs 
one h point from each diagonal, which is what parameterization /i(v) of (5.10) achieves by 
virtue of the bijection (5.11). The upshot is that 

Ap2p(u) = inf {Ap2i{h{v)) - h{v) ■ u} 
(5 24) veintw 

= inf {— /i(v) • u| = —h(u) ■ u. 

vGintW 

The last step is calculus: from the explicit formula (5.10), /i(v) • u is uniquely maximized 
at V = u. This completes the proof of Theorem 5.2. □ 

Proof of Theorem 5.3. Since \Xn\i = n and by (5.22), 
Ap2i{h)= lim n-^\ogE[e^"k=o9(Tx,w)+h-x„^ 

n— >oo 

= lim logE[e^kZo9{Tx^y^)+h{uyx„^ _ ^^^^^ + h2 = -/i2(u) + /i2. 

n— >oo 

On the other hand, by (5.23), 

log J^p(z)e5(-)+'^-+^"(-'^) = log^p(z)e^?('")+'^(")-^+^"(-'^) - /i2(u) + h2 

z z 

= -/l2(u) + /l2. □ 
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6. Limits of ratios of point-to-line partition functions 

Armed with the hmits of Theorem 4.1 and the large deviation bound of Remark 5.5, 
we prove convergence of ratios of tilted point-to-line partition functions. With the tilt 
parameter h = (/ii,/i2) G and Z^^v defined as in (4.1), let 

^i{N) = Yl e^'^""""^ Z^,v for G N and |n|i < N. 

v:\v\i=N 

This is the same as (5.3) with a general initial point. Recall definition (5.12) that associates 
a velocity u{h) = {u{h), l—u{h)) to a tilt h, and definition (3.14) that associates a parameter 
9{v) to a velocity v. 

Theorem 6.1. Fix < p < oo and let i.i.d. Gamma(p) weights {wx}xei^\ given. For 

A G (0,/)), let {^^,'n^,C^,w) be the gamma system constructed in Theorem 4-1- Then for 
h = (/ii, /12) G X € N X Z+, 2/ G Z+ X N, ¥-a.s., 

(6.1) hm /"V^ =r?^("('^)) and hm = C^"^")^ 

In other words, the limit of ratios of point-to- line partition functions tilted by h is equal 
to the limit of ratios of point-to-point partition functions in the direction u{h): 

yh ry 

hm — ; r — = hm 



•2) 



Z 7 

1 1- x,(N) ,. ^xAm,n) 

and hm — ; r — = hm 



provided m/n —?■ u{h)/{l — u{h)). We see the duality between tilt and velocity from Remark 
5.5 again. 

In Section 7 the limits of ratios from Theorems 4.1 and 6.1 give convergence of polymer 
measures to random walk in a correlated random environment. The remainder of this 
section proves Theorem 6.1. 

Proof of Theorem 6.1. We prove (6.1) for the horizontal ratios (first limit). Begin with a 
lower bound and let 5q> 



Z^ 

X 



'x,{N) _ ^^^^ Zx-e]_,v Zx,v 



^a::-ei,(Af) ^.|^|j=jv x-ei,(7V) 

Zx,v 



X] Qx-ei,(Ar){^W-|x|i + l - V) 



v:\v\^=N Zx-e^,v 



, \ ^ „/, r -lA /AT \i ^x.(m,N—m) 

> 2^ Q^_,^_(^){Xjv_|,|^+i = (m,Af-m)} 



m:\m—Nu(h)\<Ndo ' ' 
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Above we introduced a tilted quenched point-to-line polymer measure 



(6.3) QiiN)M = ^e'^<^--^yi.-y^ n < 

for paths x, from xq = y to the line |x7v-|j;|Ji = N. 

Apply construction (4.11) to the gamma system {^^ , rj^ , (■^ , w) to define partition func- 
tions and associated polymer measures with northern boundary weights {?/i*Ar_m+i}i<j<m+i 
and eastern boundary weights {Cm+i j}i<j<N-m+i- Recall the dual exit points (4.14)- 
(4.15). By an application of Lemma 9.1 (to the reversed rectangle), 



^x-ei,{m,N-m) -^^.ej ,(m+l,Af-m+l) (*ei > 0) 



lA r.* ^ n^ x,{m+l,N-m+l) 



— Qx,(m+l,N-m+l){'^ei > 0} 



,(m+l,Af-m+l)l''ei ^ ^\ 

x—ei,{m+l,N—m+l) 



Qx,{m+l,N-m+l){'^ei > ^} Vx ■ 



The last equality came from Lemma 4.3. Note the notational distinction: (fj\ is the 



tilted point-to-line polymer measure, while Q'^^y is the point-to-point polymer measure 
with boundary parameter A. 
We have the lower bound 

,,,, .Jzi""^ ^ E QU,w{^^Hx.U+i = (m,Ar-m)} 

(0-4) x-ei,{N) m:\m- Nu{h)\<N So 

^ Qx,{rn+l,N-rn+l){'^*ei > 0}^x- 

Let < A < 6'(u(/i)). Define parameter M / oo by iV(l - u{h)) = M\['i(6'(u(/i))). Let 
{m,n) = x + ([M^'i(/3 — A) J, [M\I'i(A)J), a velocity essentially characteristic for {X,p). As 
m varies in the sum on the left side of (6.4), let (mi, ni) = (m + l,A^ — m + 1). Since 
is strictly decreasing, if we fix (5o > small enough, there exists Eq > such that, for large 
enough A^, 

n - ni > M^i{X) - M'I'i(6'(u(/i))) -N6o-2> Meq 
and nil — rh > Meq. 

Following the idea of Lemma 4.4 and (4.22), 

P[Qx,i^+i,N-^+i){t:. > 0} > e-^^^°^] < P[g^,(^,^){4 > Meo} > e''^^^^] 

< g-cieoM 

Since there are 0{N) m-values, Borel-Cantelli and (6.4) give, for large enough n, 
;r > ^'(1 - e-'^'"*OQte„(iV){l^iV-|.K+i - A^u(/.)I < A5o}. 

^ ^x-ei,(Ar) 
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By the quenched LDP (5.20) for the point-to-hne measure, the last probabiUty tends to 1. 
Thus we obtain the lower bound: 

lim ^'^^^ >Vx as we let A e{u{h)). 

For the upper bound we first bound summands away from the concentration point of the 
quenched measure. 

h-{{m,N-m)-x) ry 

^x,(m,N —m) 



5^ p-hiZ^^ 
m:\m- Nu(h)\>N So x-ei,{N) 

E h-{(m,N-m)-x) V 
^ ^xAin,N —in) 


m:\m- Nu{h)\>N So 

< u7,_eiQ^,(;v){l^iV-|x|i - Nnih)\ > N6o} 0. 
For the remaining fractions we develop an upper bound. 

^x,{m,N-m) ^ ^x,(m+l,N~m+l)(^e2 ^ ^) 



Zx-ei,{m,N-m) -^^.gi ,(„i+l,7V-m-|-l) (*e2 > 0) 

1 

< 



7A 

^a;,(m+l,Af-m+l) 



Qx-ei,im+l,N-m+l)i^e2 > ^x-ei. 



{m+l,N-m+l) 



Vx 



Qx-ei,(m+l,N-m+l) 



Combining these, 



< E Qx-e.,iN)i^N-\xW+i = {m,N- m)} 

x-ei,{N) m:\m-Nu{h)\<N5o 

^x 



^ Qx-ei,{m+l,N-m+l)i^ei > 0} 

where the o(l) term tends to zero P-a.s. Proceed as for the lower bound, this time choosing 
9{u{h)) < \ < p to show that the Q'^-probability above vanishes exponentially fast. This 
completes the proof of Theorem 6.1. □ 

7. Limits of path measures 

As in Section 4, fix /> € (0, oo) and assume that i.i.d. Gamma(p) weights w = {wx ■ x G 
Z^} are given on a probability space (0,(5, P). Let Zu,v be the point-to-point partition 
functions defined in (4.1), with associated quenched polymer measures 

(7.1) Qu,v{x.} = ^ — n ^^l' x.eYiu,v 

Let point-to-line polymer measures be defined as before in (5.17) or (6.3). 
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For A G {0,p), let {S,''^ , f]'^ , , w) denote the gamma system of weights constructed in 
Theorem 4.1. In this environment, define RWRE transitions on lA by 



, and IT ' (x, X + 62 ) = 

ix+ei ' Sx+e2 ix+ei ' Sa; 



(7.2) T,^^(x,x + ei) = -j^^^^^ — and ^"''^(x, x + 62) - 

Let P"''-*' be the quenched path measure of the RWRE started at 0. It is characterized by 
the initial point and transition 

(7.3) P^'\Xo = 0) = 1, P'^'^Xk+i = y\Xk = x) = 7t^'\x, y). 

We wrote P'^'^ instead of P'^'-^ because the quenched distribution is a function of the weights 
w, through the limits (4.3) that appear on the right in (7.2). In other words, the probability 
space has not been artificially augmented with the variables that appear in definition (2.16): 
everything comes from the single i.i.d. collection w. 

Let denote partition functions defined by recipe (3.10) in gamma system (^'^, r/^, (^'^, w). 
Adapt the notation from (3.3) in the form 



(7.4) = 
Then we can rewrite transition (7.2) as 



Cx+e2' Z = 62- 



'x,x+z _ V'^O.x-- 

'^x,x+ei ~^ ^x,x+e2 (-^0,x+ei) ^ ~^ i^0,x+e2) 



(7.5) 7T-'\x,x + z)= , '-f- = V-Tr/zA -Mei,e2}. 



In other words, this RWRE is of the competition interface type defined by (2.11) in Lemma 
2.2. The next theorem shows that these walks are the limits of the polymer measures on 
long paths, both point-to-point and point-to-line. 

Theorem 7.1. The following weak limits of probability measures on the path space (Z^)^+ 
happen for F-a.e. w. 

(i) Let < X < p and suppose (m, n) — > 00 m the characteristic direction of parameters 
{X,p) as defined in (4.2). Then (5o,(m,n) converges to P"'-'^. 

(ii) Let /i G R^. Then as N ^ 00 the tilted point-to-line measure (Af) converges to 

pw,e(u(h)) 

Proof. Fix a finite path xq^m with xq = 0. Then (m,n) > xm for large enough (m,n), and 

7 M-1 M-1 A 

(7.6) 

^ ^ M-1 

= H Tr'"'\xi,Xi+i) = P'"'^{Xo,M = xo,m}. 

We applied limits (4.3) and used property Wx = Vix+ei + Cx+e2 °f gamma system 
{(^^ , T]^ , (^^ , w) from Theorem 4.1. There are countably many finite paths and these de- 
termine weak convergence on the path space. Hence P-a.s. limits (7.6) give claim (i). 
The proof of (ii) is the same with limits (6.1) instead. □ 
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Figure 3. Illustration of the proof of Theorem 7.2. The thickset RWRE path 
avoids the disk of radius SN^^^ (dark grey small disk) but enters the disk of radius 
bN^^^ (light grey large disk) centered at (kui, nn) = N\i\^p. 

The RWRE p^'^ has the fluctuation exponent of the 1+1 dimensional KPZ (Kardar- 
Parisi-Zhang) universality class: under the averaged distribution, at time n, the typical 
fluctuation away from the characteristic velocity of (A,p) is of size ii?/^. The reason is 
that the RWRE is close to a polymer and we can apply fluctuation results for the shift- 
invariant log-gamma polymer. Below E denotes expectation over the weights w. Recall the 
characteristic velocity UA,p from (3.13). 

Theorem 7.2. There exist constants Ci, C2 < oo such that for G N and b > Ci, 

(7.7) EP"''^{|Xiv - iVuA,p| > 6iV^/^} < C2b-^. 
Given e > 0, there exists 5 > such that 

(7.8) lE^ EP'^'^HXn - N\ix J < SN^/'-^} < e. 

Proof. For each let {m,n) = (^[cN"^i{p — A) J, [cA^^i(A)J) where c > is fixed large 
enough so that m A n > 2N. Define < k < 1 by = c(^'i(/9 - A) + ^'i(A)). Then up 
to errors from integer parts {nm, nn) = Nux^p. (See Figure 3.) 

Fix {m,n). We couple the RWRE p^^-^ with the polymer that obeys the quenched 
distribution Qo'(mra) defined by applying the construction (4.11) to the gamma system 
(^'*', r/^, 1^'^, w). In other words, the boundary weights ri^ and (^^ are on the north and east, 
the bulk weights come from w, and the distribution of the weights is described by (3.7), 
with w taking on the role of ^. This is the stationary log-gamma polymer to which results 
from [27] apply. 

Define the path X. G ^o^{m,n) by letting it follow the RWRE until it hits either the north 
or the east boundary of the rectangle {0, . . . , m} x {0, . . . , n}, and then follow the boundary 
to {m,n). The next calculation shows that the quenched distribution of X_ is Qg'fmn)' 
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X, G Ilg To be concrete, let < /c < m and suppose x, hits the north boundary at 

k+n—l I k+n—1 

p-^\x=x.)= n ^f^=^>^ n < 

^ k+n—l m ^ k+n—l m 



= Y\ — n ■ n ^"^ir,) ^ = ^ANE n ^^^z ■ n (^*>^ 

0,(m,n) j=0 i=k+l 0,{m,n) j=0 i=k+l 

= Qo,{m,n){^-}- 

The last equality is the definition of <5o,'^m,n)i^-}- '^^^ equality -^o,(m,n) ^ ^o,(m,n) comes 
by applying Lemma 4.3 to a telescoping product of ratio weights. 

With c large enough, the boundary does not interfere with behavior around (ktu, ku) = 
Nux^p. In (7.7)-(7.8) we can replace EP"''^{-} with EQ^~J^^^{-}. The result follows from 
Theorem 2.3 of [27], after a harmless reversal of the lattice rectangle to account for the 
difference that in [27, Thm. 2.3] the boundary weights are on the south and west. □ 

8. The log-gamma polymer random walk in random environment 

In the previous section we saw that the limits of log-gamma polymer measures are poly- 
mer RWREs with transition (2.17), where the weights come from a gamma system with 
some parameters (A,/?). In this section we identify a stationary, ergodic probability distri- 
bution for the environment process of a polymer RWRE. We expect this stationary Markov 
chain to be the limit of the environment process when its initial distribution is an appro- 
priate gamma system (Remark 8.3 below). 

The process of the environment as seen from the particle is 

Tx„^ = (Cx„+N2>'^X„+NxZ+,Cx„+Z+xN,^X„+Z2^)- 

The state space of this process is the space f^NE of weight configurations oj = (■?,??, CiO 
that satisfy NE induction, as defined in Definition 2.3 and (2.16). 

Let < a, /3 < oo and p = a + (3 + 1. Define probability distribution /i"'^ on the space 
J7ne as follows: let the variables (r/^ei , CNe2 > ) be mutually independent with marginal 
distributions 

(8.1) rjifl ~ Gamma(a), Co,i ~ Gamma(/3), ~ Gamma(/3), i,j G N. 

The remaining variables {rix,C,x,^x-ei-e-2 '■ ^ £ are then defined by north-east induction 
(2.12)-(2.13). 

A few more notational items. Ga denotes a Gamma(a) random variable and E generic 
expectation. Let P denote the distribution of the random walk on that starts at and 
has step distribution 

^(^i) = = 1 -P(e2)- 

a + p 

Let us call this the (^j^,^j^) random walk. An admissible path is denoted by xo,n = 
(xq, xi, . . . , Xn) with = and steps Zk = Xk — Xk-i G {ei, 62}- 
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The Burke property is not valid for jj!^'^ because p 7^ a + /3, so under /x"''^ the weights do 
not form a gamma system (Definition 3.1). However, the ^ weights stih turn out to have a 
tractable distribution which we record in the next proposition. 

Proposition 8.1. Under , the marginal distribution of {Cx}xei,'^ ^-^ given as follows. 
Let {hx}xei? arbitrary bounded Borel functions on Then for n G N, 



n ^-i^- 



i.2) 



x&?,: |a;|i<n 



xi,„e(z2.)- fc=o 

X JJ E/iy(G„+^+i). 

|j/|i<n: j/^{a;o,n} 



In other words, the distribution of the ^ weights is constructed as follows: run the 
(^j^) I'andom walk, put independent Gamma(a + /?) variables on the path, indepen- 
dent Gamma(a + /3 + 1) variables off the path, and average over the walks. 

Theorem 8.2. Let the environment uj have initial distribution p"''^ on the space VLj^e of 
(2.16), and let the walk Xn obey transitions (2.17). 

(a) The environment process Tx„u} is a stationary ergodic Markov chain with state space 

(b) The averaged distribution of the walk Xn is the homogeneous (^j^, ^j^) random 
walk. 

Note the contrast in the behavior of the walk X„. According to Theorem 7.2, when 
the environment has the distribution of a gamma system of weights, the averaged walk 
has fluctuations of order n^/'^. By part (b) above, when the environment has the /x"'^ 
distribution, the averaged walk is diffusive. 

Remark 8.3 (Simulations). Suppose the environment process starts from a gamma system 
with parameters (A,p), with p > 1. Simulations suggest that then Tx„uj converges to /x°'^ 
such that a + /3 = p — 1 and (^j^, ^j+;a) = UA,p, the characteristic direction (3.13) of the 
original setting. 

Under the environment distribution /i"'^, the averaged distribution of the walk Xn is the 
diffusive (^j^, ^j+;a) random walk. Simulations suggest that under its quenched distribution 
the walk localizes, with a positive fraction of overlap between two independent walks in the 
same environment. 

Remark 8.4. We can look at the environment as seen from the walk with a more general 
boundary, instead of simply the axes. Let a = {yj}j^z be a down-right path in that 
goes through 62, and ei. That is, y_i = 62, yo = 0, yi = e\ and yi - yi-i G {ei, -62}. 
Let J = {x : 3k & N : X — {k, k) G a} he the lattice strictly to the northeast of a. Weights 
assigned to this setting are such that {^x '■ x G J} are i.i.d. Gamma(p). On the path edge 
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weights have different recipes to the northwest and southeast of the origin: 

horizontal edge northwest of 0: i < 0, yi — yi^i = ei : r/y- ~ Gamma(a + 1) 

vertical edge northwest of 0: z < 0,yj — = —62 : Cyt-i ~ Gamma(/3) 

horizontal edge southeast of 0: i > l,yi — yi-i = ei : ry^,. ~ Gamma(a) 

vertical edge southeast of 0: i > — yi-i = —62 : Cyi-i Gamma(/3 + 1) 

These weights are stationary as we look at the system centered at Xn- The proof goes along 
the same lines as given below. 

Remark 8.5 (A degenerate limit and an invariant distribution as seen from a last-passage 
competition interface). The results above require p > 1. In the limit a\0, /3\0, 
the fy~^,C~^ weights blow up. We rescale so that logarithms of edge weights converge to 
exponential random variables and bulk weights vanish. Let e > 0, p = ea + e/3 + 1, and 
consider the weights H^liflwKi >Cz^''xn) wilder the distribution /x^"'^^. The independent 
weights of (8.1) now satisfy 

(8.3) ~ Gamma(p), -qf^ ~ Gamma(ea), Cg^j ~ Gamma(e/3), i,j G N. 

We can construct the weights in (8.3) as functions of uniform variables as in (4.26). Then 
the following limits as e \, can be taken pointwise: 

- e log if] 0, -e log r/J^J /j,o ~ Exp(a) , and -e log cj j ^ Jo,j ~ Exp(^) . 
The NE induction equations (2.12) converge to the equations 

(8.4) Ix — (Xj— 62 Jx—e\) and Jx — i^Jx—ei -^a;— 62)^' 

The RWRE transition probability converges to a deterministic transition: 

The limit leads to an invariant distribution for a last-passage system. Equations (8.4) 
describe inductively the increment variables 

Ix = Co.x — Go^x-ei and Jx = Go,x — Go,x-e2 

of a degenerate last-passage model with boundary weights {Ii,o,Jo,j '■ hj € N} and zero 
bulk weights. This distribution on (/nxZ+5 «^z+xn) is invariant for the environment seen 
from the location that starts at 999 = and obeys the transition 

(8.5) ^i°x+ei = '^{Ix+ei < Jx+Si} and Vr^°i+e2 = '^{Ix+e^ > Jx+ez} ■ 

Given the environment, this defines a deterministic path ip, on Z^. We recognize in (8.5) 
the jump rule of the competition interface (2.4). 

The remainder of this section is taken by the proofs. To prove stationarity of the Markov 
chain it suffices to consider the partial environment (r/^ei , CNe2 ) Cn^ ) because the other vari- 
ables of the state are functions of these. The notation here is that ij^^ei = {'niei}ieN, and 
similarly for other cases. The next lemma proves everything in Proposition 8.1 and Theorem 
8.2 except the ergodicity. 
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Lemma 8.6. Fix n € N and an admissible path xq „ with xq = 0. Fix a finite set X C "^^ 



disjoint from (x„ + Z^) U {xA,.}o<fc<n- Let {hk}kez+ o-nd {gu}uei\ collections of hounded 

B Orel functions on M_|_. Let f he a hounded Borel function on . Then 

(8.6) 

■,/3 



n-1 



fc=0 u£l 

n-1 

P{Xo 



A:=0 



Remark 8.7. Note that the independent (^Xfc) cannot go up to = n because = 
ryx„+ei + C,x„+e2 and these belong in the future of the walk. Adding the statements over 
xo,n gives the invariance of pP^'l^ and the distribution of ^. For a fixed xo,n we get the 
averaged distribution of the walk and also the statement that when the walk looks at the 
^ weights in its past, it sees Gq-i-^- variables on its path and G^+^+i-variables elsewhere. 

Lemma 8.6 is basically a consequence of size-biasing beta variables. The formulation we 
need is in the next lemma, whose proof we leave to the reader. 

Lemma 8.8. Let the gamma variables below with distinct subscripts be independent. Then 

Gn . / „ Gn \ / „ Gr- 



E 



Ga + Gp 

a 



a + P 



• ^^^) ^(g.+.+i • ^^^) KG. + G,) 
E/(G,+i)-E<7(G;3)-E/i(G,+;3). 



Proof of Lemma 8.6. We assume that the first step of the walk is ei and calculate the distri- 
bution. Introduce functions <1> to represent north-east induction (2.12)-(2.13), specifically 
to calculate the ^ weights on the vertical line x • ei = and C weights on the vertical line 
X ■ ei = 1, for X ■ 62 > 1: 

1+62) ^ei+e2+Ne2j 

= ('$i(r/ei+e2) Ce2+Ne2) '^ei+e2+Ne2 ) i '^£1+62 > ) ^2(f?ei+e2) Ce2+Ne2) '?ei+e2+Ne2 )) • 

^ %i + Ce2 ^ 



Let /iqjS') fi be bounded Borel functions of their arguments. The first equality below im- 
plements definitions. Li the second equality below apply (8.7) to the triple {Ga,Gp, Ga+p+i) 
= (ryei , Ce2) ^61+62) and note that all other variables are independent of this triple. Let 
C^+^+i denote an i.i.d. Gamma(a-|-/3 + 1) sequence. Augment temporarily the probability 
space with independent Ga+i and Gp variables that are also independent of all the other 
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variables in /2. 

^, r „ „ 

P'^iXi = ei)/lo(Co)5(?Ne2)/l(%i+Nei)/2(Cei+Ne2)/3(Cei+N2) 

,,a,/3 



X gi^liCei+ez ^".^ 1 Ce2+Ne2> '^ei+ea+Nea ) ) 

V ^ riei+ Ce2 ^ / 

X f2\iei+e2 ~T~ 1 ^2{iei+e2 ' Ce2+Ne2: '^ei+e2+Ne2 ) ) 

V %i+Ce2 ^ %i+Ce2 ^/ 



a + /3 

X ''^ [9{^l{Ga+l-, Cea+Nea, Cei+e2+Ne2)) 

X /2(G_5 , 4>2(Gq,+i , Ce2+Ne2; Cei+ea+Nea)) ] 

" E[/io(G,+/3)] n9{Gl%^,)] E^"^' ihime,)] E^^'' [/2(CNe2)] E^^^' [h{^,,.)] 



= ^ E[/.o(G,+/3)] E[5(G^:Vl)] [/l(WeJ /2(CNe2) /sl^N^)] • 

In the second last equality, inside /i and /a we simply shift by — ei. Inside /2 variable 
G/3 furnishes Ce2- Here is the key point: at this stage the Burke property applies to the 
mappings (<I>i, <I>2) because Ga+i furnishes r?ei+e2 and thereby the parameters of the input 
weights satisfy (a + 1) + /3 = p. The beta size-biasing put us back into the setting of a 
gamma system. Thus (<^i,<I>2) outputs two independent sequences. The first one denoted 
is i-i-d. Gamma(a-|-/3-|-l) and it represents the distribution of ^Ne2- The second 
one is i.i.d. Gamma(/3), which we take to be Ce2+Ne2- the last equality we can combine 
the three ^"'^-expectations because the independence is in accordance with the definition 
of 

Standard arguments generalize the product /1/2/3 so that 

E^"^^ [P'^iXi = ei) /lo(^o)5(lNe2)-^('?ei+Nei, Cei+Ne2, Cei+N^)] 

(8.9) 

= p(ei)E[/io(G,+/3)]E[g(G^^^^+,)]ii;^"'' [F(r?Mei, CNe2, ^nO]- 

for Borel functions hQ,g,F such that the expectations make sense. Reflection across the 
diagonal gives the alternative formula where the first step is 62 instead of ei , ^^62 is replaced 
by Inei, and is replaced by 

Referring to the goal (8.6), let Zq = X \ (xi -|- Z^) and take g{^,) = rTuGXo 
can combine the ei and 62 cases into this statement, which is (8.6) for n = 1. 

^M-./S = Xi) /lo(lo) • n ■ E{Vxi+Ne„ Cxi+Ne2> L^+N^)] 

(8.10) 
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To obtain (8.6), do induction on the length n of the path. Let X' = X n (xi + Z^). In 
(8.10) take 



n— 1 n— 1 



i=l fc=l uGX'— 

Assuming (8.6) holds for paths of length n — 1, the right-hand side of (8.10) turns into the 
right-hand side of (8.6). □ 

The ergodicity claim of Theorem 8.2 is in the next lemma. 

Lemma 8.9. With initial distribution fi"'^ , the stationary process Sn = (f?x„+Nen Cx„+Ne2) 
?X„+N2) is ergodic. 

Proof. Denote a generic state by 5 = (??Nei, CNe2; ■^N^)- It suffices to show that, for any 
function / G (//"'''), the averages 

n-1 

n-i^i^^[/(5fc)] 

k=0 

converge to a constant in L^[^"'^) [26, p. 91-95]. By approximation in L^{fj,"'^), it 
suffices to prove this for a local function /, that is, a function of the variables s = 
(^4,0) Co,j) Ci,j)i,je[M] for arbitrary but fixed M G N. Let s = ip{S) denote the pro- 
jection mapping, and let the projection of the stationary process Sn be s„ = (/?(5„) = 
iVXn+ihO) ' Cx„+{o,j) , ^x„+ii,j) )i,je[M] ■ 

Process s„ is also a stationary Markov chain, with state space and invariant 

distribution v = /i"'^ o ip~^. Under v coordinates of s are independent with distributions 
ijifi ~ Gamma(a), C,oj ~ Gamma(/3) and ^ij ~ Gamma(p). 

Given the state s = (vifi, Co, j,^i,j)i,j^[M]j we compute the variables {rjxXx '■ x € [M]^} 
via north-east induction (2.12). The transition from state s to a new state goes by two 
steps: (i) randomly shift s by ei or 62; (ii) add fresh variables to the north or east to replace 
the variables lost from south or west in the shift of the M x M square. 

Precisely speaking, from s = {rjifiXo,jyCi,j)i,je[M] the process jumps to either t' or t", 
according to the following two cases. 

(a) The shift is ei and t' = (??j+i,0) Ci,j) '^i+i,i)i.je[Af] where the new independently 
chosen variables are rjM+ifi ~ Gamma(a) and Cm+ij ~ Gamma(/3) for j £ [M]. 

(b) The shift is 62 and t" = (??i,i, Coj+i; ^ij-i-i)i,je[M] where the new independently 
chosen variables are Co,A/+i ~ Gamma(/3) and ^i,M+i ~ Gamma(/3) for i £ [M]. 

The probabilities of the two alternatives are 

7r(s,tO= and 7r(s,0= • 

Vei + Ce2 Vei + Ce2 

Let 7r(s, dt) denote the transition probability of the Markov chain s„: the shift followed by 
the random choice of new coordinates to complete the square [M] x [M]. The task is to 
check that s„ is an ergodic process. 
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Two general observations about checking the ergodicity of a Markov transition P with 
invariant distribution v. (i) Suppose has a density with respect to a background measure 
A. Then it is enough to check that, for z^-a.e. x, P{x, dy) has a density p{x, y) with respect 
to \{dy) such that y) > for A-a.e. y. For then, if A is a v-a.s. invariant measurable 
set such that v^A*^) > 0, taking x ^ A'^ la 



1a{x) = P{x,A)= / p{x,y)\{dy) 
J A 



shows that A(^) = and thereby y{A) = 0. (ii) It is enough to check the ergodicity of 
some power P™. 

We show that for m = 2M + 1, 7r™(s,(it) has a Lebesgue almost everywhere positive 
density on M^m+m^^ Let S be a Borel subset of R^^+^^'. Write 

(8-11) TxS = (?/a;+{i,o)>Cx+(Oj)'Cx+{i,i) : hj G W]) 

for the shifted configuration in the M x M square. 

7r"(s,i?) = Yl i^^"' [1b(7;s)Po"{^™ = x} \ ^{Sq) = s] 

xdT?. ■.\x\\=ra 

(8 12) 

x(iT?^:\x\\=m 

On the first line above E^^ ' represents the choices of fresh coordinates while the shifts 
are in the quenched probability Po'jXm = x}. After that we conditioned on the cr-algebra 
(Figure 4) 

(8.13) Tlx = cr{r]nei, CNea; Cx, {Cij : i < X ■ ei - I OT j < X ■ 62 - I}} ■ 




Figure 4. The cr-algebra Hx- The dark black sites in the interior and the thickset 
lines on the axes denote the 77, (} variables that generate "Hx- The gray lines de- 
note {77, C} variables computed via north-east induction from information contained 
in Hx- Finally the lighter gray sites denote <^ variables independent of Hx- 
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Xm = X implies |x|i = m, and then m = 2M + 1 guarantees that Hx is large enough to 
contain the event ^{Sq) = s. The quenched probability PfflXm = x} is also ?^x-nieasurable. 
Of the variables that make up T^s in (8.11), the ^x+(i,j)^^ independent of Tij., but the 
^x-i-(i,o)'s and Cx+(Oj)'s depend on Ti^ through the equations 

(8.14) Vx+(i,o) = ix+(i,o) ^"^"^['7^^ ' « = 1, • • • , ^ 

and 

(8.15) CxMo,) = Uio,j)- '^-^^^ , j = l,...,M. 

Vx+{0J-1) + ^x+{-l,j) 

The situations for {??x-i-(i,o)} ^'^d {Cx+{o,j)} symmetric, so let us look at equation (8.14) 
closely. Tix contains variables {Cx ; rjx+ii,-!) '■ ^ ^ i^]} because these can be computed by 
north-east induction from the variables listed in (8.13), so these are taken as given in 
(8.14). Variables {(,x+{i,o) ■ ^ ^ i^]} are picked i.i.d. Gamma(p), independently of Tlx, 
while variables {Cx+{i-i,o) -1 = 2,..., M} are calculated along the way from the equations 

(8.16) Cx+{i-lfl) = Cx+{i-lfi) 2 0) ^ i = 2,...,M. 

Vx+ii-l-l) + <5X+(i-2,0) 

Regarding {(x ', ?7j.+(j _i) : i G [M]} as given parameters, equations (8.14) and (8.16) show 
that the vectors fy = (?7a;+(i,o) • ^ € [M]) and ^ = {Cx+{i,o) '■ i ^ i^]) ™ (0) oo)*^ are bijective 
functions of each other, and these functions are rational functions with positive coefficients. 
(The coefficients themselves are functions of {(^x i]x+{i.-i) '■ ^ ^ [-^]}-) Thus the Jacobians 
of these functions cannot vanish on {0,oo)^ . Consequently, from the everywhere positive 
density of (product of Gamma(p) distributions) we get an everywhere positive density /i 
for ??, for every given value of {Cx ; Vx+{i-i) ■ i G [M]}. 

This argument can be repeated to get an everywhere positive density /2 for the vector 
C = iCx+{o,j) '■ i £ for every given value of the variables specified by the conditioning 

on Tix- 

Let / denote the (everywhere positive) density of the vector {S,x+(ij) ■ h3 ^ [-^^D- With 
this notation we can write 

E^'''''[lB{TxS)\nx]= I lB{u,V,w)hiu)f2{v)f{w)dudvdw 

where the right-hand side is not a constant, but the densities fi and /2 depend also on the 
variables specified by the conditioning on l-Lx- On the right the densities are multiplied due 
to independence that comes from dependence on disjoint sets of variables. This formula 
can be substituted into (8.12) to conclude that 7r™'(s, •) has an a.e. positive density on 

(0,00)2^+*^'. □ 



9. Appendix: Auxiliary results 



This section contains a comparison lemma for partition functions, a large deviation bound 
for the log-gamma polymer, and an ergodic theorem for correctors. 
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9.1. Comparison lemma for partition functions. Let arbitrary weights {V^j^^g^a^ be 

given and define partition functions as in (2.5). For a subset A C Yl^^y, define tlie restricted 
partition function (unnormalized polymer measure) by 

\v-u\i 

(9.1) Z^,,{A) =Y. U 

Recall the definitions of the exit points (3.8)-(3.9). The restriction A = {t^i > 0} means 
that the first step of the path is ei. In other words, ZQ^^ite-i > 0) = V^Y^-^ei.xj defined for 
a; • ei > 1. 

Lemma 9.1. For m>2 and n > 1 we have this comparison of partition functions: 

Zo,{m~l,n){'tei ^ ^) ^ ^ .^0,{m-l,n) (^£2 > 0) 

'^0,(m,n)(*ei > 0) ^^(1,1), (m,n) .^0,(m,ra) (^62 > 0) 

Proof. Consider the ratio weights for these partition functions: 

_ Zo,x-ei{tei > 0) _ Zei,x-ei ~ _ Z{l,l),x-ei 
Zo^x{tei > 0) .^61,2: Z{l,l),x 

^ _ Zo,x-e2{tei > 0) _ .^ei,x-e2 ^ _ ^(l,l),x-e2 
Zo,x{tei > 0) Zej^^x Z{^i^i)^x 

On the boundary of the lattice these ratios satisfy 

Ci,j = ^ij = Ci,j and rji^i = Vi^i ^ < = for i,j>2. 

NE induction (2.12) preserves these inequalities, and gives the first inequality of (9.2). The 
second comes analogously. □ 

9.2. Large deviation bound for the log-gamma polymer. Let < a < p and let 

(C;^)C) be a gamma system of weights with parameters (a,/?) according to Definition 3.1. 
Let Zq,^ be the partition function defined by (3.10) in this gamma system, with the corre- 
sponding point-to-point quenched polymer measure 

, *cxit Hi 

(9.3) ^3o,.K} = ^(n-r.LM)( n ^^1)' --^no,. 

Let the scaling parameter > 1 be real valued. Let (m, n) € denote the endpoint of 
the path. Measure the deviation from characteristic velocity by 

(9.4) KN = \m- N'^i{p - a) I V I n - iV^'i(a) |. 

Lemma 9.2. Let hn be defined by (9.4). Let 6 > 0. Then there are constants < 6i,c,ci < 
oo such that the following estimate holds. For {m,n) G , > 1 andu > {IV ckn^/ SN), 

(9.5) ^[Qo,im,n){tei > > e"'^" ] < e-^^^ 

Same bound holds for te2- The same constants work for {a,p) that satisfy < a < p and 
vary in a compact set. 
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Proof. Let /3 < a and take two gamma systems: ??°,C°) with parameters (a,p) and 
i^iV^iC^) with parameters {f3,p). Couple them so that they share the ^-variables, and 
r]i < 7]^ and (x > C hold. This can be achieved by imposing these same conditions on the 
variables in part (c) of Definition 3.1, and then noting that the inequalities are preserved 
by (3.f). Let Z°' and be partition functions computed in these two systems. 



(9-6) 

^ 0,{m,n) I J Vifi 



exit 1 m+n 

^ - - -.ndr)( I 



*cxit 2^ s , m+n 



'■exit 



n 



0,{m,n) j=i 'i,0 

In the bounds below, X = X — KX denotes a centered random variable. Recall the mean 
(3.11). Let 5i > 0. From (9.6) 

[u] 



F[Qo,(m,n){ie, > u} > e"^^^] < p| log ryf „ - logr^-^ ) > d,u\ 

i=i ' 

^^■^^ +^{ log<(m,n)-log^oV,n) > (N -m)(^I/o(a) -^'0(/3)) 

+ n(^o(/> -13)- ^'o(p - a)) - 2-51^}. 
Standard large deviations apply to log- gamma variables, so 3c2 > such that 

L«J >, 

j;(logr?fo-logr?-o) ><5in <e-^^^ 
1=1 ^ 



(9.8) 



Taylor expand to second order the ^'o-differences inside the last probability in (9.7). Keep- 
ing 6 > fixed, pick 5i > and a — /3 > small enough and c < oo large enough. Then 
for another small constant C3 > 0, the probability simplifies to 



(9-9) H log <(m.n) - log ^0"(m,n) > ^3^} < e^^". 

The bound comes again from i.i.d. large deviations, by virtue of (3.12). □ 

9.3. Ergodic theorem for correctors. With a bit of extra effort and with future use 
in mind we prove this ergodic theorem more generally than required for this paper. Fix a 
dimension d G N. Let {Q, (3,P) be a probability space equipped with a semigroup {Tx+y = 
Tx o Ty) of measurable maps : — )• $7 for x G Z^J.. Generic points of J7 are denoted by 
bj. Assume P invariant and ergodic under {Tx).j.^-^d : that is, FoT^^ = P, and if T^^A = A 

Vx G Z'^ then F{A) G {0, 1}. Let TZ = {a : i = 1, . . . ,d} denote the set of admissible steps, 
which are precisely the standard basis vectors of W^. Admissible paths {xk)^^Q satisfy 
Xk - Xk-i G TZ. 

Let FifixT^— T-Mbea corrector, by which we mean these properties: 
(i) Vz G 7^: F(w, z) G ^^(P) and EF(a;, z) = 0. 
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(ii) The closed-loop (or cocyle) property: if {a^^j^^Q and {x'^}^^q are two admissible 
paths such that xq = x'q and Xn = x'^, then 



n—l n— 1 

/ ^ 

fc=0 i=0 



Define the path integral of F by 

f{x, u)) = Y^ F{Tr,^uj, Xi+i - Xi), {x, uj) eZ'lx n, 

2 = 

where (xj)"^Q is any admissible path from xq = to Xn = x. /(0,w) = 0. The closed-loop 
property ensures that / is well defined. 

We make the following additional assumption. 

(9.10) Iim Iim max - V |F(r^+j^a;, z)[ = Mz&Tl. 

(5\0n-i>oOa;g2d .1^1 <„ n 

+ 0<j<n<3 

A similar assumption was useful in [24, 25] in studies of polymers. If, for each z TZ, the 
variables {F{TizUJ, z)}i^z+ are i.i.d., then by Lemma A. 4 of [25] a sufficient condition for 

(9.10) is 

(9.11) 3p >d :mF{uj,z)\P] < oo. 

Our application of Theorem 9.3 is to the corrector F^ in (5.14). By Corollary 5.1 this 
satisfies the i.i.d. condition and even has an exponential moment. Thus hypothesis (9.10) 
is satisfied by F^ in (5.14). 

Theorem 9.3. Let F be a corrector and satisfy (9.10). Then 

(9.12) lim max = 

■>''->-°° xeZi:\x\i=n ri 



As an auxiliary result towards the main theorem, we prove a limit for averages over 
rectangles of any dimension. The following result is a discrete version of Lemma 6.1 of [18]. 

Theorem 9.4. Let F be a corrector. Let r £ [d], < ai < hi for I <i < r, and let 

(9.13) Bn,r = {x G : [nffljj < Xj < [nbij for I < i < r, Xr+i = • • • = = 0}. 
Then 

(9.14) lim — y = ¥-a.s. 

It is enough to consider the case = 0, for the general case is obtained by successive 
differences and sums of such cases. Then to simplify notation we take 6j = 1. We separate 
a part of the proof as a lemma. 
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Lemma 9.5. Take Ui = and hi = 1 in (9.13). Let I < j < r < d and g : [0,1]'" 
continuous. Then 

(9.15) lim ^ V g{n~\xu...,Xr))F{T^oj,ej) = P-a.s. 

n^oo 77,' ^ — ^ 



Proof of Lemma 9.5. Fix j. The pointwise ergodic theorem [19, Thm. 6.2.8] gives 

hm 4 F(r,u;,e,) = /i(^^). 

We wish to show that h is invariant under each shift Tg . . By the closed-loop property (now 
for j G {1, . . . ,r}) 



^ n— 1 n— 1 
^ ■■■ -^(^fciei+fc2e2+---+fcrer^) Cj) 



fci=0 fcr=0 

n— 1 n— 1 n— 1 n— 1 



"I ^ ■ ■ ■ X/ ^ ••• -^(^fciei+fc2e2-|---|-ri.e,-|---|-fc^er'^; 6») 



n 

fcl=0 fcj_i=0 fcj + l=0 fcr=0 



^ n—l n—l n—l n—i 

^ X] ■ ■ ■ X] ■ ■ ■ X^ ^(^fciei+A:2e2+-+0-ej + ---+A:rer'^, Cj) 

A:i=0 fcj_i=0 fcj+i=0 A;r=0 

^ n—l n—l 
+ — X^ ■ ■ ■ X] -^(^fciei+fc2e2+-+fcrer(7'eiW),ej). 



A:i=0 kr=0 

The closed loop above is undertaken for each fixed r^^eiH i-A:j_iej_i-i-fcj+iej+iH hferer^- The 

two paths are {ej, 2ej, . . . , nej,nej + Cj} and {ei,ei + ej, Cj -|- 2ej, . . . , Cj -|- ncj}. 

The first sum converges to /i(w), the last one to h{Te^u). By the pointwise ergodic theo- 
rem the first sum on the right converges to because it has only n*""^ terms. Consequently 
all terms converge a.s. The second sum on the left must also vanish in the limit because it 
converges to zero in probability. We get h{u}) = h(Te-u), and conclude by ergodicity that 
h = 0. Then (9.15) follows by a Riemann sum-type approximation. □ 

Proof of Theorem 9.4- This goes by induction on r. For r = 1, rearrange: 

n ^-^ n n ^-^ n ^-^ n ^-^ \ n I 

k=0 k=l i=0 k=o ^ ^ 

An application of (9.15) gives the conclusion (9.14) for r = 1. 
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Suppose that (9.14) holds for some rG{l,...,d— 1}. Let us show it for r + 1. 

x£Bn r+1 x£B„.r fe=0 



xGB„,r k=0 

X^Bn.r x£Bn,r k = 

rf ^-^ n 71^+^ ^-^ \ n I 

x£Bn,r X&Bn,r + l 

As n — )• oo on the last line, the first sum goes to zero by the induction hypothesis and 
the second sum by (9.15). □ 

Proof of Theorem 9.3. It is enough to prove 

(9.16) lim min Zi^ili^^ > Q. 

n— >-oo xGZ[j_:|a;|i=n 

This statement applied to both F and —F verifies (9.12). 

Let (5 > and = k6/{4d) for k € For k = {ki, . . . , k^) G Z^^., define rectangles 
Bn,k = {x & : [najfc-J < x-i < [ncfc^+ij, i G [d]}. For each x € 1j\ such that |x|i = n 
we can pick a rectangle B^^x = ^(x) such that every point y G B^^x can be reached from 
X with an admissible path of at most n5 steps. Our strategy is to replace f{x,uj) by an 
average of / over this rectangle. Note that there is a fixed finite set K of vectors k such 
that the above choices can be made from {i?„^k : k € for all n and all = n. 

For every x such that = n and every y G Bn^x fix a path from x to y such that the 
steps 61,62, ... ,6^ are taken in order. Then for any such pair x,y, with designated path 



f{x,uj) = f{y,uj) - y F{Tx^uj,Xi+i 



m—l 

Xi 
i=0 

>/(y,w)-V| max y \F{Tu+izUJ,z)\\. 

z£TZ +11— 0<i<nd 

The error term is independent of x,y. Average over y G -Bn,x and then take minimum over 
= n: 

. f{x,uj) ^ . 1 f{y,^) 

mm 



> min 



\x\i=n n kgX \Bn,]i\ ^ n 



y\ max - y \F{Tu+izU},z)\\. 

zGTZ +11— 0<r<nd 
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As n ^ oo, the first term on the right vanishes by Theorem 9.4. After that let (5 — > and 
assumption (9.10) takes care of the last term. Bound (9.16) has been verified and Theorem 
9.3 thereby proved. □ 
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